Extremal extensions of entanglement witnesses : Unearthing new bound entangled 

states 
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In this paper, we discuss extremal extensions of entanglement witnesses based on Choi's map. The 
constructions are based on a generalization of the Choi map due to Osaka, from which we construct 
entanglement witnesses. These extremal extensions are powerful in terms of their capacity to detect 
entanglement of positive under partial transpose (PPT) entangled states and lead to unearthing 
of entanglement of new PPT states. We also use the Cholesky-like decomposition to construct 
entangled states which are revealed by these extremal entanglement witnesses. 
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I. INTRODUCTION 

Quantum entanglement plays a central role in quan- 
tum theory from a conceptual as well as a practical point 
of view. On the conceptual front, entanglement is inti- 
mately connected with the notions of non-locality and 
violation of Bell's inequalities, which lie at the heart of 
the way the quantum mechanical description of the world 
differs from the classical one. On the practical front, 
quantum entanglement is essential in providing a com- 
putational advantage to quantum computers over their 
classical counterparts [1]. 

The quantum state of a physical system with a finite- 
dimensional complex Hilbert space H (which is cither 
pure or mixed), is represented by p £ B{H) which is a 
positive semi-definite hcrmitian operator with unit trace. 
The set of states forms a convex set and the extremal 
points of this set are pure states which are operators of 
rank 1. 

For composite quantum systems, the Hilbert space is 
the tensor product of Hilbert spaces of the individual 
systems. Consider a bipartite quantum system with its 
state space given by Ha <8> Hb where Ha and Hb are the 
Hilbert spaces of individual quantum system. A bipartite 
state pab € B(Ha <£> Hb) is said to be separable if it is 
possible to decompose it as follows: 



PAB 



j'=i 



pb 3 ; Vj V] > 0, 



i; (i) 



where pAj and pBj are states in systems A and B respec- 
tively A bipartite state p S B(Ha <B> Hb) is said to be 
entangled if it is not separable i.e. it cannot be expanded 
in the form given in Equation (1). 

For pure quantum states, the characterization of states 
as separable or entangled, is easily achieved by comput- 
ing the entropy of the reduced density operator of the 
subsystems. However, such a characterization for mixed 
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states is a non-trivial problem. While a large volume of 
work has appeared on this issue in the past two decades, 
the problem still remains open. In the absence of a "final 
solution" , explorations into finding new classes of mixed 
entangled states and new ways of constructing them, is 
useful and provides insights into the classification prob- 
lem. 

An important method to detect the entanglement of 
quantum states is to construct entanglement witnesses, 
which arc positive maps that are not completely positive 
[2]. Consider a positive map Ha defined on the Hilbert 
space Ha of system A. If this map is not completely 
positive, then for some Hb the map h,A®I acting on Ha® 
Hb will not be positive. Therefore, there will be a state 
Pab for which (Jia <8> I) Pab will be negative. However, 
this cannot happen for a separable state defined in (1). 
Therefore, such a pab has to be entangled. 

Partial transposition was the first such entanglement 
witness which was used to unearth the entanglement of 
pure as well as mixed states [2, 3]. While negativity 
under partial transpose (NPT) indicates entanglement, 
positivity under partial transpose is both necessary and 
sufficient only for 2 (g) 2 and 2 <g> 3 systems. The proof re- 
lied on the earlier works of Woronowicz [4], Arvison [5, 6] 
and St0rmer [7, 8] where they show that in dimension 2, 
all positive but not completely positive maps are decom- 
posable (can be written as a combination of a completely 
positive map and a transposed completely positive map) . 
A corollary of their result shows that, for dimensions 2(g) 2 
and 2 <S> 3, a state is separable if and only if it remains 
positive under partial transpose. However in higher di- 
mensions there could be maps which are positive but not 
completely positive and which may not have a connection 
with partial transpose. This means that there could be 
states which are positive under partial transpose (PPT) 
and are still entangled. Such states are called PPT en- 
tangled states. Obviously for such states, there exists a 
positive but not completely positive map revealing their 
entanglement [4, 9-12]. 

Choi provided the first example of an indecomposable 
map which later led to the detection of entangled states 
beyond partial transpose [13]. A number of other discrete 
examples are also available [14-22]. However there is no 
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systematic way of characterizing such maps beyond two 
dimensions. A few methods have been recently suggested 
for generating examples of positive maps which are not 
completely positive [23, 24]. However, there is no straight 
forward procedure to verify if they are indecomposable 
or not. Generalizations of Choi's method have also been 
adopted to produce indecomposable maps which in prin- 
ciple have the potential to reveal entanglement of new 
quantum states [25]. 

The present paper is an attempt in this direction, 
where we have constructed extremal extensions of known 
entanglement witnesses in order to unearth new classes of 
bound entangled states. We have constructed examples 
of PPT entangled states which are revealed by extremal 
Choi type maps considered by Choi and Lam [26], and 
later by Osaka[16], and have exploited the Cholesky de- 
composition in a nontrivial way. The Cholesky decom- 
position has been used by Chrusciriski ct al. [27] in the 
context of PPT states; however our analysis goes much 
beyond that. Moreover, given any extremal positive but 
not completely positive map, we can use our method to 
generate new classes of extremal maps. In particular, 
for the 3 (g> 3 bipartite system, we have constructed PPT 
states whose entanglement is revealed by our map but is 
not revealed by Choi's map. For the positive maps which 
are not completely positive, the question of extremality 
has been settled only for a few maps and the problem 
of determining the structure of such maps and their ex- 
tremal or non-extremal nature is in general difficult. In 
this regard, we believe that our generalization adds new 
insights into the class of extremal entangled witnesses. 
Further, for a given map which is positive but not com- 
pletely positive, it is not always trivial to find the PPT 
entangled states revealed by the map. In our case we 
were able to find a family of such states. 

The material in this paper is arranged as follows: In 
section II we discuss the bi-quadratic forms and their 
connection with positive maps, where the notions of ex- 
tremality and dccomposability arc discussed as impor- 
tant ingredients for the purpose of classification of en- 
tangled states. In section III we explore the possibility 
of extremal extensions of Choi's map. We present Os- 
aka's map as an extremal extension of Choi's map and 
go on to construct other extremal extensions. We then 
turn to constructing examples of bound entangled states 
whose entanglement is revealed by Osaka's map and by 
our new extremal extension of Choi's map. At the end 
of this section we show that the bound entangled states 
that we have constructed are robust and can take a cer- 
tain amount of noise before they lose their entanglement. 
Section IV offers some concluding remarks. 



II. BI-QUADRATIC FORMS AND EXTREMAL 
MAPS 

Given a non-negative polynomial P(x\,--- ,x n ) > 
0, V(xi, • ■ • ,x n ) £ M. n of degree d, the question whether 



P can always be written as a sum of squares of poly- 
nomials has been around for a long time. Minkowski 
conjectured that in general the answer should be 'no'. It 
was proved by Hilbert that, except for three exceptional 
cases (n = 1, d arbitrary; n arbitrary, d = 2; and one 
non-trivial case n = 2, d = 4), there always exist pos- 
itive semi-definite polynomials which cannot be written 
as a sum of squares of polynomials. However, his proof 
was by an indirect method and did not provide actual 
examples of such polynomials. For a survey and devel- 
opment of the problem, see Rudin [28]. The generalized 
version of this problem on rational polynomials, is known 
as Hilbcrt's 17th problem, and for which the answer is 
'yes'. 

The first counterexample was constructed by Choi [13] 
where he considered a positive semi-definite bi-quadratic 
form F„(X : Y) (each term having degree four), with 
six variables, divided into two sets denoted by X = 
{xi,x 2 ,x 3 }, and Y = {2/1,2/2,2/3} given by: 

F^X:Y) = {xlyl+xlyl + xlyl) 

-2(xix 2 yiy2 + x 2 x 3 y 2 y 3 + 2:3^12/32/1) 
+fi(xjy 2 2 + x\y\ + x\y{) (2) 

where \x > 1. 

Choi proved that for \i > 1, this bi-quadratic form is 
non-negative definite but cannot be written as a sum of 
squares of quadratic forms [13]. Later it was shown by 
Choi and Lam [26] that it is also true for fj, = 1. Choi's 
method has been modified and extended, and differ- 
ent examples of such positive semi-definite bi-quadratic 
forms were found. Among these, the results by Os- 
aka [16], Cho ct. al. [14], and Ha [15, 17, 29] are impor- 
tant. Later generalizations of these methods for gener- 
ating such forms in arbitrary dimensions were developed 
by Chruscihski and Kossakowski [25]. 



A. Connection with positive maps 

The intimate connection between a positive map and 
positive semi-definite bi-quadratic forms was also discov- 
ered by Choi [13]. Before describing the connection we 
provide a few definitions. 

Definition : A hcrmiticity preserving map h : 
M n (C) — ► M n (C), is said to be a positive map, if it 
maps positive semi- definite operators to positive semi- 
definite operators. Here M„(C) denotes the set of all 
?i x n complex matrices. 

Definition : A positive map h is said to be fc-positive 
if the extended map 

h <g> l ft : Mft(C) <g> M„(C) — ► Mft(C) ® M n (C), 

is positive, we here 1ft denotes the identity mapping on 
the auxiliary space Mft(C). The map h is said to be 
completely positive if the above extensions are positive 
for all k > 1. 



3 



The connection between the maps and bi-quadratic 
forms can be established as follows. Consider a hermitic- 
ity preserving linear map 



S 



(3) 



We can construct the corresponding bi-quadratic form 
F(X : Y) as 



F(X :Y) = (Y\S{X ■ X T )\Y) 



(4) 



where X = (xi ■ ■ ■ x m ) and Y = (yi • • • y n ) ■ 

On the other hand, let F(X : Y) be a bi-quadratic 
form. Notice that, it is a quadratic form with respect 
to Y (as well as X). So we can write it in the form 
(Y|Ax|y). Thus we get a map which takes any one- 
dimensional projection X.X T to Ax- Using linearity 
and hermiticity, we can extend it to a map which pre- 
serves hermiticity. It was shown by Choi that, given any 
positive semi-definite form, the corresponding map is a 
positive map and vice- versa [13, 30]. 

There is thus a bijective relation between the set of 
positive semi-definite forms and positive maps between 
matrix algebras. The property of complete positivity can 
also be translated easily. If a map is completely positive, 
the corresponding bi-quadratic form can be written as a 
sum of squares of quadratic forms and vice versa. Put 
differently, if a map is positive but not completely pos- 
itive, the corresponding bi-quadratic form will be pos- 
itive semi-definite but can not be written as a sum of 
squares of quadratic forms. Thus each such form gives 
rise to a unique map between the space of real sym- 
metric operators, which can be trivially extended to the 
set of all hermitian operators, and then to all operators. 
This also connects with the work of Arvison [5, 6] and 
St0rmer [7, 8] who were exploring the set of positive maps 
between C*-algebras. Since then, other examples and 
classes of such maps have been discovered. 

By the above correspondence, the Choi quadratic form 
given in Equation (2) leads to the following map for 3x3 
matrices. 



*c(M) 




an + M«33 — ai2 -oi3 
—02i a 22 + ^an ~a 23 
— 031 -°32 a 33 + fj,a 2 2. 



, (5) 



with fi > 1. From the quadratic form (2), exchanging 
the X and Y variables, we can get another map, which 
is given by 



«n ai2 a 13 ' 
$q (p) : I 021 a 2 2 a 2 3 I 

(G31 a 32 0-33/ 



an 



r M a 22 
a 2 i a 22 
a 3 i 



a\ 2 

f- va 3 3 

032 



033 




(6) 



Our interest in these positive but not completely pos- 
itive maps is because of their ability to detect entangle- 
ment of quantum states. For the maps that are to be 
used as entanglement witnesses, two notions, namely de- 
composability and extremality are very important. We 
define these notions below. 

Definition : A positive but not completely positive map 
is called decomposable, if it can be written as a sum of a 
completely positive and a completely co-positive map. 

This property was first discussed by Woronowicz [4]. 
Since decomposable maps are obtained by combining a 
completely positive map with a transposed completely 
positive map, it is clear that they are weaker than partial 
transpose in terms of their ability to detect entanglement 
and therefore are not of interest. The interesting point 
however is that, given a map which is positive but not 
completely positive, there is no standard way to check if 
it is decomposable or not! 

Since the set of positive maps is a convex set it can be 
described by its extremal elements. Therefore, it is most 
natural to study extremal positive maps. From the point 
of view of detecting entanglement, any extremal map is 
more powerful than the maps which are internal points 
of the set of positive maps [16, 26]. Choi and Lam define 
an extremal map using the corresponding bi-quadratic 
form [26] as follows. 

Definition : A positive semi-definite bi-quadratic form 
F is said to be extremal if, for any decomposition of F = 
F\ + F2 where F^s are positive semi-definite bi-quadratic 
forms, Fi = XiF, where Ai, A 2 are non-negative real 
numbers with Ai + A2 = 1. 

It was shown by Choi and Lam that in the case [i = 1 , 
the form Fi defined in Equation (2) is extremal. 

Since the set of positive semi-definite forms is a convex 
set, it is enough to identify the set of such extremal forms. 
At this stage it is useful to change the notation to the 
original Choi-Lam notation for the bi-quadratic forms 

f x\ x 2 £ 3 N 



and we therefore denote F(X : Y) as F 



Vi 2/2 2/3 



III. EXTREMAL MAPS AND BOUND 
ENTANGLED STATES 

Osaka extended the result of Choi and Lam, and gen- 
erated a class of extremal maps [16]. Osaka's map 
$0(^,2/1 z ) is defined as 

(an 012 ai 3 
021 022 023 I 
031 032 033 



an + xa 33 

-0,21 

-a 3 i 



a 22 



-ai 2 
+ yau 
-032 



033 



-013 

a 23 

+ za 22/ 



,(7) 



where x,y, z > and xyz = 1. Osaka showed that this 
class of maps is extremal [16]. 

Generalizing beyond Osaka construction, we define a 
class of extremal bi-quadratic forms as follows: 
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Let F = F { Xl X2 2:3 ) be an extremal positive semi- 
\Vi 2/2 2/3/ 

definite bi-quadratic form. 

For a set of non-zero positive real numbers a,b,c we 

define 



G 



Xl x 2 x 3 

yi 2/2 2/3 



F 



ax 1 bx2 cx 3 

2/1 V2 2/3 



(8) 



/ £pi X 3 \ 

It turns out that the form G is positive 

\yx V2 yz) 

semi-definite and extremal. 

Proof : We first prove the positivity. Let us assume 
that the proposition is not true and there exists real num- 

(Pi P2 Pz\ 
9l 92 93 J 

But then by definition F { Pl P2 P3 ) < 0; for a set of 

\9i ° 2 #V 

real numbers p'i,p 2 ,p 3 , 9i> 92, 93 where p{ = api, p' 2 = 
bp 2 and p' 3 = cp 3 . This contradicts the assumption that 
F is a positive semi-definite form for all real values of 
Xi's and yj's. 

For extremality, let us assume G = G\ +G 2 , where Gi 
and G 2 are positive semi-definite bi-quadratic forms. 
Notice that 



Xl x 2 x 3 
yi 2/2 2/3 



G 



Gi 



£1 


■('2 










a 


b 










m 


2/2 


2/3 J 










X2 




/ Xl 






u 


b 




+ G 2 ( « 


6 






2/2 


2/3 J 


\2/i 


.'72 


2/3 j 



as we have assumed. Define two positive semi-definite 
forms 



Xl x 2 x 3 
2/1 2/2 2/3 



Xi X2 X3 

a b c 

2/1 2/2 2/3 



for i = 1,2. We can now write, 
F = Fi + 



F 2 



However, F is extremal. Therefore, Fj = A;F, i = 1, 2, 
and Ai and X 2 are positive real numbers with Ai +A2 = 1. 

Xl x 2 x 3 N 



Xl £2. £3. 

Thus G, ( a 6 c 
,2/1 2/2 2/3 



Xl x 2 x 3 
2/1 2/2 2/3 



: A 4 F 

= G, 

= A,F 
= AiG. 



2/i 2/2 2/3 



Hence 



a £l 5^2. c £3 

a c 

2/1 2/2 2/3 

axi bx 2 cx 3 
2/i 2/2 2/3 



Since Ai, A2 > and Ai + A2 = 1; the form G is an 
extremal form. 

The maps corresponding to the bi-quadratic form G 
defined above are positive maps which are extremal. 
This construction is extendable to higher dimensions 
without any further work. This means that given 
any extremal positive semi-definite bi-quadratic form, 



F = F 



Xl x 2 



,2/i 2/2 • 
01, a 2 , ■ ■ ■ ,a n 'i G 



for any non zero positive real 



F 



is posi- 



aiXi a 2 x 2 ■ ■ ■ a n x n 

2/1 2/2 ••• Vn 
tive semi-definite and extremal. 

We now turn to the map corresponding to the extremal 
bi-quadratic form G from Equation (2). After working 
out the details, the map turns out to be 



Xu Xi 2 X13 

3>(a, b,c) : | x 2 i x 22 x 23 
,X 3 i x 32 x 33l 



' a 2 xn 



c 2 x 33 



-abx 2 i 
—acx 3 i 



—abxi 2 —acxi 3 
b 2 x 22 +a 2 xn —bcx 23 |(9) 
-bcx 32 c 2 x 33 + b 2 x 22 . 



where a, b 1 c 7^ 0. The map $(a, b, c) is an extremal posi- 
tive map which may not be completely positive. For the 
values of a, b and c for which it is not completely positive, 
it can act as an entanglement witness. 

We now construct a set of PPT entangled states for 
which the above map acts as an entanglement witness. 
Consider a density operator for a 3 ® 3 system defined by 
two parameters t and x. 



/ l+t 
t 
i 



p{x,t) 



-At 




\ 






x 



x 





x 









x 







x \ 


a; 




x 
a; 



1 
\ 
0l/ 
(10) 

This p is a unit trace density operator for t > and 
< x< 1. 

The action of the map $(a, b, c) on the density operator 
p(x,t) leads to the transformed density operator p x \ 



p'{x,t) = (<f>(a,b,c)®l 3 )p(x,t). 



(11) 



We compute eigenvalues of p(x, t)' and use the negativ- 
ity of the least eigenvalue as an indicator of entanglement 
of p(x, t). It is useful to note that the map $(a, b, c) with, 
a = b = c = 1 reduces to Choi's map $^(1) wnne f° r 
other values of a, b and c it is still an extremal map with 
a potential to reveal entanglement of quantum states. 

A computation of eigenvalues of p'(x,t) reveals that 
for this example, the maps $(a, b, c) have more potential 
than Choi's map in unearthing the entanglement of PPT 
quantum states. The results are displayed graphically 
in Figures 1 and 2. We take the parameter values to 
bea=l + |,fe = c= l and calculate the minimum 
eigenvalues of p(x, t) for the range x & [0,1] and t € (0,1]. 
The results are displayed in Figure 1 . The curved surface 
denotes the minimum eigenvalue of p{x, t) after the action 
$(1 + |,1,1). To show the power of this map clearly, we 
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display a section of the above graph where we fix the 
parameter x = ^ and plot the minimum eigenvalue as a 
function of t. We compare our result with Choi's maps. 
The result is shown in Figure 2. The continuous line 
here denotes the minimum eigenvalue corresponding to 
the action $(1 + = , 1, 1) while the other two dashed lines 
are the minimum eigenvalues corresponding to & c and 
It turns out that approximately after x = 0.604428, 
the minimum eigenvalue becomes negative under $(1 + 
|, 1, 1) while the minimum eigenvalues under $^(1) and 
$p 7 (l) still remain positive. It is only after x crosses the 
value 0.66 that Choi's maps begin to detect entanglement 
for this class of states. Therefore, for p(x, there is 
a clear window of x values where the entanglement is 
revealed by<&(l + |,l,l) and is not revealed by any of 
the Choi's maps. 

The map <&(1 + |, 1, 1) was chosen as a representative 
example. In fact the class of maps $(a, b, c) can reveal the 
entanglement of a large class of PPT entangled states and 
therefore provide a genuine extremal extension of Choi's 
maps. 



0.60 




FIG. 1. The plot of least eigenvalues of p'(x,t) as a function 
of x and t. The curved surface corresponds to the case where 
p'(x,t) was generated by the action of $(1 + |,1,1) upon 
p(x,t). The middle plane is the plane xy — 0, given a re- 
ferral plane for highlighting the negativity of the eigenvalues 
represented by the curved surface. 



A. PPT entangled states detected by Osaka's map 

We now address the question of constructing PPT en- 
tangled states for which Osaka's map acts as an entan- 
glement witness. Although the three-parameter family of 
maps due to Osaka as described in Equation 7 has been 
defined, and is known to be positive but not completely 
positive, there has not been an explicit construction of 
PPT entangled states whose entanglement is revealed by 
this class of maps. 

We set up a computer search over the PPT entan- 
gled states and employ the Cholesky decomposition [31] 
to selectively scan the states in the Hilbert space. We 
also make an intelligent use of Choi- Jamiolkowski isomor- 
phism developed in [9, 13, 32, 33] to check at each stage of 




FIG. 2. The section corresponding to t = ^ of Figure 1 
is displayed here. The blue line corresponds to the map 
<E>(1 + |, 1, 1) while the other two curves correspond to the 
two Choi's maps. The window of values of x (approximately 
x > 0.6025 and x < 0.66) where the map $(1 + |, 1, 1) is able 
to reveal the entanglement of p(x, t) and where Choi's maps 
do not reveal the entanglement is clearly visible. 

our search that the state remains entangled. This method 
of searching for entangled states for a given positive but 
not completely positive map is in fact more general and 
can be tried for other maps too. 

According to the Cholesky decomposition, every den- 
sity matrix p of a quantum system can be decomposed as 
p = TtT, where T is an upper triangular matrix [31]. If 
p is strictly positive (i.e. all eigenvalues are greater than 
zero), then the corresponding Cholesky Decomposition is 
unique. For simplicity, we restrict ourselves to the case 
of real states where elements of the density matrix are 
real and in this case the Cholesky decomposition reduces 
to p = T f T, where t denotes the transpose operation. 

The Choi Jamiolkowski isomorphism developed in [9, 
13, 32, 33] provides a simple one-way test of entanglement 
for a given positive map which is not completely posi- 
tive. Consider a composite system Hi <8> H2 where both 
the subsystems arc of dimension d. Using the standard 
basis |i) in both the H\ and H2 we define a maximally 
entangled state 

1 d-l 
* i— 

Given a positive map $ that is not completely positive 
we define an operator 

d-l 

= — 53 |*)<<| (g. ^(Ii><i|). (13) 

* i— 

Given a density matrix p defined on H\®H2 the operator 
W provides a sufficient condition for entanglement 

Tr(lFp) < p Entangled (14) 

This is a one-way condition and Ti(Wp) > does not 
imply that the state p is separable. The condition for 
entanglement given in (14) helps us in quickly identifying 
states whose entanglement is revealed by the map $ and 
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wc employ this condition in our search for PPT entangled 
states revealed by Osaka family of maps. 

We first construct a upper triangular matrix T. We 
further restrict to those states which are invariant under 
partial transpose by imposing the condition (T'T) PT = 
T'T (where PT denotes the transpose with respect to the 
second system) . This amounts to generating a non-trivial 
solution of the equation 

(T t T) PT - T*T = 0. (15) 

It is not always possible to find a non-trivial solution to 
the above equation. However, if we begin with a sparse 
matrix T we can hope to find a solution to the above 
equation. In our search we also at every stage impose 
the following condition 

Tr(W$ o p)<0 m 
Tr{W$ c ip) > 

This means that we restrict ourselves to those PPT states 
whose entanglement is revealed by Osaka's map but is not 
revealed by Choi's map. The above methodology guides 



us in our computer search process to look for classes of 
states whose entanglement is revealed by the Osaka fam- 
ily of maps. 

Using this method and employing a computer search 
protocol, we construct an example of a PPT entangled 
state for a3®3 system. The upper triangular matrix T 
given by 

/ Vwy ooo^o oo \ 

y 2 + 5?/ 

3y 1 + 2/0 

0000 00 

0000^0 i+» 
vio viu 

00 l + 2y l + y 

0000 00 

0000 00 

\ 0/ 

(17) 

leads to a one parameter family of density operators, pa- 
rameterized by a positive parameter y. 



10y 2 













2/(2 + 5y) 











32/(1 + 2/) 





y 2 





2/(2 


+ 5y) 




































32/(1 + 27) 











2/(2 + by) 





(2- 


h5y) 2 

















y(2 + 5y) 













2 + 6y + by 2 











(l + y) (l + 2y) 



















(l + 2y) 2 





(l + y)(l + 2y) 











3y(i + y) 













{i + y? 

























(l+»)(l + 2y) 





(1 + 2/) 2 





3y(i + y) 













(l + y)(l + 2y) 











(i + y) 2 . 



where N = 10 + 36y + 57y 2 is the normalization fac- 
tor such that Tr(p(y)) = 1. We apply a one parameter 
subfamily of Osaka's map defined in (7) $o (X^, j) to 
the family of states p{y) and compute the eigenvalues of 
the resultant operator. We do a similar computation of 
the eigenvalues of the operator which is obtained by the 
action of Choi's maps and <&qI for comparison. 

In Figure 3 the least eigenvalue is plotted as a function 
x and y. Here the curved surface denotes the minimum 
eigenvalue corresponding to the $o (l, x, The middle 
plane denotes the xy plane which is placed to indicate 
the place when the surface becomes negative. In Fig- 
ure 4 we have taken a fixed value of x. The eigenvalues 
are plotted along the vertical axis and y varies along the 
horizontal axis. We apply the map $o (l, 6, g) to this 
state p{y) and plot the minimum eigenvalue which is de- 



noted by the continuous curve in Figure 4. The dashed 
curve denotes the minimum eigenvalue achieved by the 
Choi's map. The plot highlights that approximately af- 
ter point y = 0.326402 the minimum eigenvalue under 
<&o (1,6, becomes negative. Thus there is a range of 
values, where Osaka's map can identify more PPT en- 
tangled states while the Choi's map fails to do so. 



B. Robustness analysis 

We now consider the robustness of the states given in 
Equations (18) and (10). Let p be an arbitrary entangled 
state. We consider the convex combination of p with a 
maximally mixed state. For this case we consider the 
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FIG. 3. Application of Osaka's map <E>o (l,£, 5) to the den- 
sity operator p(y). The two independent variables are x and 
y, and the vertical axis denotes the eigenvalues of p(y) under 
the map. The curved surface represents the variation of the 
minimum eigenvalue. The plane in the center is the plane 
xy — 0, which highlights the portion of the surface with neg- 
ative eigenvalue. 




0.369284 



FIG. 4. Application of Osaka's map with x = 6. The hor- 
izontal axis is y and eigenvalues are along the vertical axis. 
The continuous line denotes the variation of eigenvalues un- 
der $0 (l, 6, |). The dashed line corresponds to Choi's map. 
It clearly shows that approximately after y = 0.326402 on- 
ward, PPT entanglement is revealed by $o((l, 6, ^). How- 
ever Choi's map can reveal entanglement approximately from 
y — 0.369284. This shows the superiority of Osaka's map over 
Choi's map for this instance of p(y). 

following convex combination; 

p'{e) = £ -l Q + {l-e)p; 

and try to detect the range of e for which p' is entangled. 
I9 denotes the identity matrix of dimension 9. Typically, 
the map which detects entanglement of p is used on p' as 
well. 

We begin with p(y) of the example (18). Using the 
process previously described, we obtain the new state 

p'(e,y) = |Iq + (1 - e)p(y). 

For y = I, p{y) is an entangled state, whose entan- 
glement is revealed by $0. We use the map $0 on the 



family of states p' (e, |) and can see that there is a con- 
tinuous range of e for which p' {e 1 1) remains entangled. 
The change in eigenvalues is shown in Figure 5. It shows 
that approximately up to e = 0.047, the state remains 




FIG. 5. Application of Osaka's map $0 on the convex com- 
bination p' (e, |) = |l9 + (1 — e)p (|). e is plotted along 
the horizontal axis, and eigenvalues along the vertical axis. 
The line shows the minimum eigenvalue as a function of the 
robustness parameter e. The state remains entangled approx- 
imately up to e = 0.045. 

We now use the same procedure for p(x, t) of the ex- 
ample in (10). The family of states is given by; 

p(e, x, t) = -I 9 + (1 - e)p(x, t). 

We use P (jq,Jq), which is an entangled state and its en- 
tanglement is revealed by 3> (l + A, 1, l). Now the iam- 
ily p' (e, jq, ^) is dependent on e. We plot the minimum 
eigenvalue of this family as a function of the robustness 
parameter e. Figure 6 shows that up to approximately 
e = 0.012 the state remains entangled. 



0.0m 

0.003- 




FIG. 6. Application of $ (I + 1, l) on p (e, ^, ^). e is 
plotted along the horizontal axis, and eigenvalues along the 
vertical axis. The line shows the minimum eigenvalue as a 
function of the robustness parameter e. The state remains 
entangled upto approximately e = 0.012. 
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IV. CONCLUSIONS 

This work is an exploration in the context of finding 
bound entangled states whose entanglement is revealed 
by witnesses based on positive maps that are not com- 
pletely positive. We have managed to extend the con- 
struction to a class of bound entangled states using the 
Cholcsky decomposition whose entanglement is revealed 
by Osaka's map acting witness. Furthermore, we 



have generated a family of extremal extensions of Choi's 
original map and shown that these extremal extensions 
are capable of revealing the entanglement of new classes 
of entangled states. We are extending this work to de- 
velop a more general framework where extremal move- 
ment in the map space is tracked down to a similar move- 
ment in the space of bound entangled states. Those re- 
sults will be presented elsewhere. 
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